Trilinear Higgs couplings in the two Higgs doublet model with CP violation 
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We carry out a detailed analysis of the general two Higgs doublet model with CP violation. We 
describe two different parametrizations of this model, and then study the Higgs boson masses and 
the trilinear Higgs couplings for these two parametrizations. Within a rather general model, we 
find that the trilinear Higgs couplings have a significant dependence on the details of the model, 
even when the lightest Higgs boson mass is taken to be a fixed parameter. We include radiative 
corrections in the one-loop effective potential approximation in our analysis of the Higgs boson 
masses and the Higgs trilinear couplings. The one-loop corrections to the trilinear couplings of the 
two Higgs doublet model also depend significantly on the details of the model, and can be rather 
large. We study quantitatively the trilinear Higgs couplings, and show that these couplings are 
typically several times larger than the corresponding Standard Model trilinear Higgs coupling in 
some regions of the parameter space. We also briefiy discuss the decoupling limit of the two Higgs 
doublet model. 

PACS numbers: 12.60.Fr, 14.80.Cp, ll.30.Er 
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I. INTRODUCTION 



The Higgs niechanism [ij of spontaneous electroweak 
symmetry breaking is a necessary ingredient of the Stan- 
dard Model (SM), which is crucial for its internal con- 
sistency. The search for the Higgs boson is, thus, one 
of the major tasks for the experiments at the upcom- 
ing Large Hadron Collider (LHC). In order to confirm 
the Higgs mechanism as the origin of spontaneous break- 
ing of SU{2)l X U{1)y gauge symmetry in the SM, not 
only must the Higgs boson be discovered, but also its tri- 
linear {X^^H ) and quartic {X^^hh ) ^^^^ couplings must 
be measured in order to completely reconstruct the Higgs 
potential. Furthermore, one must also be able to mea- 
sure the couplings of the Higgs boson to gauge bosons 
and fermions. 

Here we shall be concerned mainly with the trilinear 
self couplings of the Higgs boson. In the Standard Model, 
there is only one trilinear self coupling of the Higgs boson 
which can be written simply in terms of the Higgs boson 



mass Mh as 
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where v = 2Mw / g — 246 GeV is the vacuum ex- 
pectation value of the neutral component of the Higgs 
doublet, Mw is the mass of the , and g is the 
SU{2)l gauge coupling, respectively. Several extensions 
of the SM, such as the minimal supersymmetric standard 
model (MSSM) |1, and the general two Higgs doublet 
model (2HDM) [31 , have a more complicated Higgs struc- 
ture. In these models there are several trilinear Higgs 
couplings, having more complicated dependence on the 
underlying masses. It is a challenging task to measure [3| 
these trilinear couplings at the LHC. On the other hand a 
linear collider could possibly offer much better prospects 
of measuring these trilinear Higgs couplings [5, 6, 7, 8, 9]. 

At the proposed International Linear Collider (ILC), 
if the Higgs boson is not too heavy, the trilinear Higgs 
coupling can be measured via the Higgs boson pair pro- 
duction e+e^ W^*TiW~*v HHvv. A precise mea- 
surement of the trilinear Higgs self coupling will also 
make it possible to test extended Higgs models, which 
have a different structure of the Higgs potential, and 
hence different trilinear Higgs couplings, as compared 
to the SM. The 2HDM is the simplest, yet a very gen- 
eral, model with an extended Higgs structure, which 
leads to various distinct physical effects. In particular, 
the model can easily accommodate additional CP vio- 
lation [13, [m [H, beyond what is generated by 
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the Kobayashi-Maskawa (KM) mechanism ■ Further- 
more, since there is a distinct possibihty of measuring 
precisely the Higgs boson self couplings at the ILC, there 
is a motivation to study the radiative corrections to the 
trilinear self couplings of the Higgs boson. 

On the other hand, there are a number of parameters 
in the potential of the two Higgs doublet model, includ- 
ing those associated with CP violation, which determine 
the masses and CP properties of the model. In spite 
of the complicated nature of the Higgs potential, it has 
been shown in the 2HDM with CP conservation that one- 
loop corrections to the lightest CP-even Higgs boson self 
coupling are in general of a non decoupling nature, and 
can give rise to 0(100%) deviations from the SM predic- 
tion [1^. This happens even when all other couplings of 
the lightest Higgs boson to gauge bosons and fcrmions 
are in good agreement with the SM prediction. 

In this paper we shall study in detail the trilinear Higgs 
couplings in the two Higgs Doublet Model (2HDM) with 
CP violation. We shall assume that the underlying gauge 
group is the SM gauge group. After spontaneous break- 
ing of the SM gauge symmetry, the Higgs spectrum of 
the model consists of three neutral Higgs bosons and two 
charged Higgs bosons. In the CP-conserving version of 
the model, two neutral Higgs bosons (/i", H^) are CP 
even, whereas one neutral Higgs boson {A'^) is CP odd. 
In the CP-conserving case, there are, thus, six allowed 
trilinear Higgs couplings which can be labelled as Xhhh, 
>^hhH, >^hHH, >^HHH involving the CP-even Higgs bosons, 
and XhAA, Xhaa, all even in the number of A'^, involving 
the CP-odd Higgs boson. When CP is not conserved, 
the neutral Higgs bosons do not have a definite CP, and 
there is no such constraint on the couplings involving the 
A" Higgs boson. There, are, thus, a total of ten trilin- 
ear Higgs couplings in the case of the 2HDM with CP 
violation. 

Since this is one of the simplest models which goes be- 
yond the KM mechanism of CP violation, it is important 
to study the impact of CP violation on the Higgs self 
couplings in this model. We shall, therefore, consider 
the 2HDM with explicit CP violation in the Higgs sector 
in this paper. This explicit CP violation is introduced 
through appropriate complex parameters in the poten- 
tial of the 2HDM. Furthermore, we shall assume that 
this explicit CP violation cannot be transformed away 
by a redefinition of the Higgs fields. Thus, one of the ob- 
jectives of the present study is to determine the effects of 
this explicit CP violation in the Higgs sector on various 
Higgs self couplings in the two Higgs doublet model. In 
our study of the model, we shall find it convenient to keep 
the masses of the two lightest neutral Higgs bosons fixed, 
and then study the dependence on other parameters of 
the model, which will determine the amount of explicit 
CP violation in various Higgs self couplings, and also de- 
termine the mass of the heaviest Higgs boson. In this 
way, we do not lay emphasis on the heavy Higgs sector of 
the theory, whereas at the same time we can exhibit the 
wide range of values that the trilinear Higgs couplings 



can assume. In particular, we wish to emphasize that, 
in contrast to the MSSM, the heavy mass effects in a 
general 2HDM do not decouple. In order to demonstrate 
this we shall calculate the one-loop corrected Higgs boson 
self-couplings using the method of effective potential in 
a general 2HDM with explicit CP violation as described 
above. This is in contrast to Ref. jlBI, where all couplings 
and mass parameters were assumed to be real, thereby 
precluding the phenomena of explicit CP violation in the 
Higgs sector as considered in this paper. 

The plan of this paper is as follows. In Section |TT] we 
describe the most general two Higgs doublet model with 
CP violation. We then discuss the spectrum of the Higgs 
bosons of the model, and the constraints on the param- 
eters when there is CP violation in the model. Here we 
also discuss how only certain sets of parameters of the 
Higgs potential can lead to physically consistent models, 
and discuss different ways in which these parameters can 
be specified. We delineate the regions of the parameter 
space where there is no CP violation, corresponding to 
which one of the neutral Higgs bosons is odd under P. In 
Section Uni we derive the tree level trilinear couplings be- 
tween the neutral as well as the charged Higgs bosons of 
the model, and discuss the correspondence of these cou- 
plings with the trilinear Higgs couplings of the minimal 
supersymmetric standard model. Here we also study nu- 
merically the trilinear Higgs couplings, and discuss the 
domain of the parameter space that is compatible with 
various theoretical and experimental constraints. 

In Section HVl we use the method of one-loop effective 
potential to calculate, as a first step, one-loop corrections 
to the Higgs boson masses. In Section |Vl we then calcu- 
late the one-loop corrections to the trilinear self couplings 
of the Higgs bosons. In Section lVTl we carry out a detailed 
numerical study of the one-loop corrected trilinear Higgs 
couplings, discuss the magnitudes of different contribu- 
tions, the dependence on the scale, as well as the decou- 
pling limit, and show that these couplings can be several 
times larger than the corresponding Standard Model tri- 
linear Higgs coupling in some regions of the parameter 
space. We summarize our results and conclusions in Sec- 
tion |Vlll Some of the analytical calculations used in our 
analysis are described in appendices. 

II. THE TWO-HIGGS-DOUBLET MODEL 

The general 2HDM with the underlying gauge group 
SU{2) L X U{1)y is obtained by extending the Higgs sector 
of the SM with a second SU (2) l Higgs doublet with weak 
hypercharge Y = 1. Thus, the model contains 4 complex 
scalar fields which are arranged as SU{2)l doublets as 
follows: 

$,= (^^0 (Y^+l), z = l,2. (2.1) 

Using these two Higgs doublets, the most general renor- 
malizable potential for the 2HDM which is invariant un- 
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der the SU{2)l x U{\)y gauge group can be written as more appropriate to define the weak states 



A4($I$2)(*^$l) + 2 [A5(*I$2)' + ll.C. 
{ [A6($I$i) + A7($^$2)] («'I$2) + h.C.} 



(2.2) 



where A^ (i = 1, . . . , 7) are dimensionless parameters and 
the subscript "tree" denotes that (j2.2|) is a tree-level po- 
tential. We note that A^ (z = 1, . . . , 4) are real, whereas 
Xi (i = 5, . . . , 7) are in general complex parameters. Sim- 
ilarly, rriii and rn|2 real, whereas is in general 
complex. We note that the terms proportional to Ae and 
A7 have to be constrained, since th is p otential does not 
satisfy natural flavor conservation [l6|, even when each 
doublet is coupled only to up-type or only to down-type 
quarks. 

When the neutral components of the two Higgs dou- 
blets $1,2 acquire vacuum expectation values (VEVs), 
ipi2 = vi,2/V^, the gauge group SU{2)l x U{1)y breaks 
down to C/(l)cm, whereby three of the eight real fields 
in (|2.ip are absorbed by three of the four gauge bosons 
of SU{2)l X U{1)y, which become massive in the pro- 
cess, leaving behind a massless photon. We can then 
parametrize the two Higgs doublet fields in (|2.ip as 



-^{vi + + ■iX^) 



1, 2, 



(2.3) 



where we have chosen the VEVs of the neutral Higgs 
fields to be real, and absorbed the relative phase be- 
tween the two VEVs in the parameters TO22 ? A5 , Ag and A7 
of Vtree- When we substitute the parametrisation (|2.3p 
of the Higgs fields in the potential (|2.2p . there will be 
cubic terms in the Higgs fields arising from the quartic 
couplings Xi, which will give rise to trilinear couplings 
among the Higgs fields. It is these trilinear couplings 
that we shall study in detail in this paper. 

As discussed in the Introduction, in the case of 2HDM 
with CP conservation, after spontaneous breakdown of 
the gauge symmetry, we are left with two CP-even Higgs 
bosons h^,H^, a CP-odd Higgs boson , and a pair of 
charged Higgs bosons H^. However, with CP violation 
the neutral Higgs bosons h^,H^,A^ mix, and it is then 



(2.4) 



to describe the neutral Higgs sector of the CP violating 
two Higgs doublet model. In (|2.4p we have defined 



7/3 = - sin /3 xi + cos (3 X2, 
G" = cos/?xi + sin/?X2, 



(2.5) 
(2.6) 



where is the would-be neutral Goldstone boson, and 
tan (3 — V2/V1 is the ratio of the vacuum expectation 
values of the two Higgs fields. We shall take tan/3 as 
an independent parameter. The charged Higgs fields and 
the corresponding charged Goldstone boson are likewise 
defined as 



sin /3(pf + cos P (ff , 



G ^ cos (3 tpi -fsin/?c/?2 



(2.7) 
(2.8) 



respectively. 



A. Basis rotation 

The tree-level mass squared matrix of the neutral Higgs 
bosons can now be defined as 



M 



(2.9) 



where, after differentiation, all fields are set equal to zero: 
Vi — V2 — V3 = = G° = G^ = 0. The physical neu- 
tral Higgs states Hi , which are the eigenstates of the mass 
squared matrix (|2.9p . are then obtained by a rotation R 

H = Rri, ri = R^H, (2.10) 

or, more explicitly 

Hi = RijTjj, rjj = RijHi. (2.11) 

The rotation matrix R which diagonalizes (|2.9p . 

RM^R^ = A^Lg = diag(M2, M|, Af|), (2.12) 

with Ml < M2 < M3, can be parametrized as 



R = i?3 i?2 ^1 



'10 

cos a3 sin as 
, — sin Q!3 cos as 



cos a2 sin a2 

1 
— sin a2 cos a2 



cosai sinai 0> 
— sinai cosai 
1, 
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Cl C2 
-(Cl S2 S3 + Si C3) 
-Cl S2 C3 + Si S3 



Si C2 S2 
Cl C3 - Si S2 S3 C2 S3 
-(Ci S3 + Si S2 C3) C2C3 



(2.13) 



where Ci = cosai, Si — sincti. Note that Eq. (|2.12p can 
be inverted as 

{M%j = Rk^MlRkJ. (2.14) 

k 

By symmetry, = (A^^)jj, and, in the general case, 

this matrix is seen to contain 6 independent parameters. 
These may be taken as the three masses, and the angles 
of the rotation matrix. 

As discussed in Ref. , there are three limits in which 
there is no CP violation in the 2HDM, corresponding to 
which one of the neutral Higgs bosons is odd under P. 
These limits can be characterised in terms of the angles 
a2 and as of the rotation matrix (j2.13p as follows: 

Hi is odd: ^2 ±7r/2, (2.15) 
H2 is odd: 03 = ±7r/2, (2.16) 
i73 is odd: a2 = 0, a, = 0. (2.17) 

We note that in the latter limit, we have ai — a + tt/2, 
where a is the conventional mixing angle in the CP- 
even sector of the MSSM Any CP violation in the 
2HDM will also depend on the Yukawa couplings. In the 
2HDM II (Model II), where the down-type quarks cou- 
ple to $1 and the up- type quarks couple to $2, various 
measures of "maximal" CP violation both in the Higgs- 
vector boson and Higgs-quark sectors are discussed in 
p^ . These maxima occur in the "bulk" of the a2-Q;3 
space, typically for |q!2| = 0(7r/4) and la,] = 0{n/4:). 

B. Neutral Higgs boson masses 

Minimizing the tree-level potential according to 




and using the resulting minimization conditions to elim- 
inate mf^ and 77132 ' ^"^^ obtains the elements (|2.9p of the 
tree-level mass squared matrix 



Mil 


— v\\i + vlv 


+ ^Re {3vf As - 
27;i 


vl Ay), 


MI2 


= vl X2 + 


+ ^Re{-vl Xe H 


h 3vl Ay) 


Ml, 


= v^Re [-A5 + Tr—{vl Ag 4 

2viV2 


^vlXr)], 


MI2 


= viV2[Re (A3 - 
+ f Re (vl Ae + 


f A4 + As) - H 

vl Ay), 





Ml, = -i-wlm [7;2 As + 2t;i Ae], 



■^23 = -5«Im[wi As + 27;2 Ay], (2.19) 
where we have defined 

vi=vcosf3, V2^vsmf3, /3 e (o, |) , (2.20) 

and 

7.= -^ Re (777^2), (2.21) 

2viV2 

with t;^ — vf + vl = (246 GeV)^. The squared masses 
Mf of the physical neutral Higgs bosons are obtained as 
the eigenvalues of the mass squared matrix (|2.9p . These 
eigenvalues are solutions of a cubic equation which in- 
volves the parameters A^. However, only some set of val- 
ues of the parameters A.^ lead to consistent solutions. In 
order to identify physically consistent models, we shall, 
following the approach of Refs. [H,[23,[2l|, specify phys- 
ical (tree- level) masses, instead of the Ai, as parameters. 
Depending on how much constrained a model we want 
to consider, we shall do this in two ways, which we shall 
denote "approach (A)" and "approach (B)": 

(A) In this case, which corresponds to Ae = Ay = 0, 
two elements oi are related via tan/3 and we 
cannot therefore take all three masses as indepen- 
dent. Instead, we take the two lightest neutral 
Higgs boson masses, together with tan/3 and the 
three angles (ai, 02,03) defining the mixing matrix 
R of Eq. (|2.10p . as independent parameters. The 
third (heaviest) Higgs boson mass can then readily 
be determined. In general, the elements i?i3 and 
i?23 of the rotation matrix R must be non-zero in 
order to have CP violation. For consistency, the 
derived quantities Im As and Im 777^2 must be non- 
zero. 

(B) Here we take Xq, Ay 7^ 0. In this case, we take all 
three neutral Higgs boson masses, tan /3, together 
with the mixing matrix R of Eq. (I2.10|) . ImAs, 
Re Ag and Re Ay as the input parameters. 

In either case, specifying the charged Higgs boson mass, 
Mh±, as well as the bilinear parameter /x^ = v'^v, we 
can determine the remaining Aj through a set of linear 
relations. They are, thus, unique. Details are given in 
Appendix 1X1 

We note that in general, there exist multiple minima in 
the 2HDM potential [l^jllll- However, with our choice of 
input parameters, including Higgs squared masses, and 
these being positive, the minimum we are working in is 
a global one and hence stable [23j . 
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III. TRILINEAR SELF COUPLINGS OF HIGGS 
BOSONS AT THE TREE LEVEL 

In this Section we shall define the trilinear Higgs self 
couplings, and obtain explicit expressions for them in 
terms of the parameters of the tree-level Higgs potential 
of the 2HDM. We shall then study the behavior of these 
couplings as functions of the various parameters of the 
model. 



A. Trilinear Couplings of Neutral Higgs Bosons 

The trilinear self couplings of the neutral Higgs bosons 
are defined as 



mmdiTk' 



(3.1) 



which are most easily obtained from the corresponding 
derivatives of V in (|2.2p with respect to the weak fields rji. 
The derivatives with respect to Hi in (|3.ip are obtained 
by noting the useful relation 



drij d 



d 



which can be used to go from the rji basis to the physical 
Hi basis. 

When there is CP violation, the trilinear couplings 
among the physical Higgs bosons will involve elements 
of the rotation matrix R j2^ . [25j . We can then write 
the trilinear couplings in terms of the derivatives of the 
potential (|2.2p with respect to rji and the elements of the 
rotation matrix R as 



m<n<o— 1,2,3 



drjnidrfndrjo 



E 



Ri'mRj'nRk' o 



(3.3) 



m<n<o— 1,2,3 



where the indices m, n, o refer to the weak field basis, and 
the * denotes a sum over permutations P, {i',j',k'} — 
P{i, j, k}, which gives rise to a factor of nl for n identical 
fields. We now proceed to obtain these couplings in an 
explicit form. 

At the tree level, when expressed in terms of A^, the 
derivatives in Eq. (|3.3p are rather simple. The trilinear 
couplings arnno amoug the weak fields rji can be written 
as (in units of —iv) [24. [2^: 

iiii = |(cos/3 Ai + sin/?Re Ag), 

aii2 = I (sin /3 Re A345 + 3cos/3Re Ae), 

aii3 = — i [cos /3 sin /3 Im A5 + (1 + 2cos^ /3) Im Ae], 

ai22 = 5 (cos /3 Re A345 + 3sin/3Re A7), 

ai23 = — ImAg — cos/3 sin/3(ImA6 + ImXr), 



0-133 = 5{cos/?(sin^ Ai + cos^ /?Re A345 — 2Re A5) 

+ sin/3[(sin2/3- 2cos^ /3) Rc Ae + cos^ /3Re A7]}, 
0222 = 5(sin/3A2 +cosf3ReXr), 

0-223 — —^[cosP sin/3Im A5 + (cos^ /? + 3sin^ /3)Im A7], 

^233 = 5 {sin /3(cos^ /? A2 + sin^ (3 Re A345 — 2Re A5 ) 

+ cos /3 [sin^ /3 Re Ae + (cos^ /3 - 2 sin^ /3) Re A7] } , 

0333 = i (cos (3 sin P Im A5 — sin^ /3 Im Ag — cos^ /3 Im A7 ) , 

(3.4) 



where 



A 



345 



A3 + A4 + A5. 



(3.5) 



In order to elucidate the compact notation in Eq. 
we explicitly write the trilinear couplings Am and A112 
in Appendix [Bl 

Rather than studying all the trilinear couplings of the 
Higgs bosons, we shall here focus on the two couplings 
Am and A112 involving the lightest Higgs boson Hi. In 
the special case when Ag = A7 = 0, and if in addition 
0^2 = 0^3 = 0, these couplings take the simple form 



Am = -3iv{cl cpXi + 3/3X2 



(3.6a) 



+ cisi(ci Sf3 + si C/3)Re A345}, 
A112 = w{3ciSi(cic^Ai - S1S/3A2) 

- [cfsp - sfcp + 2ciSi(ciC/3 - SiS/3)]ReA345}. 

(3.6b) 

In this limit of Ag = A7 = and 0:2 = as = 0, we note 
the following: 

• For small values of tan/3, both Am and A112 are 
determined by Ai and ReA345. This is because for 
small tan/3, terms containing sp = sin/3 vanish. In 
actual practice, due to the constraints imposed by 
i3-physics, the small tan /3 limit is not reached. 

• For large values of tan/3, both Am and A112 are 
determined by A2 and ReA345. In this case terms 
containing c/3 = cos /3 vanish. 

Here, the combination Re A345 is related to via 

Mfj± = + ^v^ih - Re A345 - Re (Ag + A7)]. (3.7) 

It is instructive to compare the couplings (|3.6p with 
the Standard Model trilinear coupling (jl.ip . While the 
SM trilinear coupling is given in terms of only one pa- 
rameter, the mass Mh of the Higgs boson, those of the 
2HDM depend on several parameters, and may actually 
pass through zero, even for the special case studied in 
Eq. (|3.6p . In the general case, we find it convenient to 
study the dimensionless ratios of the couplings 



A 



111 



\SM ■ 



6^ 



Ai 



12 



\SM ■ 



(3.8) 



where for the reference SM coupling we use the mass 
Ml of the lightest Higgs boson of the 2HDM. These di- 
mensionless ratios of the couplings will be calculated and 
discussed in the following. 
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1. 2HDM-MSSM correspondence 



It is also useful to compare the trilinear Higgs cou- 
plings in the 2HDM with the corresponding ones in the 
minimal supersymmetric standard model. In the MSSM, 
where there is no CP violation in the Higgs sector at the 
tree level, there are six trilinear Higgs couplin gs. The 
couplings corresponding to Am and A112 are [gI l26| 



>^hhh = - „ , , /I cos(2a) sm(Q! + p), 



\hhH 



2 cos 9w 
igMz 



2 cos 6w 
cos(2q!) cos(a + /3)] 



2 sin(2Q;) sin(Q; + /?) 



(3.9a) 



(3.9b) 



where Mz is the mass of the Z, Ow is the Weinberg 
angle, and g the SU{2)l gauge couphng. The angle a 
is the mixing angle in the CP even Higgs sector of the 
MSSM. In the limit a2 — ^ and 0, it corresponds 

to ai - 7r/2 of the 2HDM. 

In the MSSM the trilinear Higgs couplings are con- 
trolled by the SU{2)l gauge couphng g. In contrast, 
the trilinear Higgs couplings in the 2HDM arise from the 
quartic terms of the Higgs potential when the SM gauge 
symmetry is broken. It is easy to see that in the above 
limit the couplings in the two models are simply related 
by the following correspondence 



Ai — A2 
Re A345 



1(9 



2 



4V.y i 9 



+ 9% 



together with 



As = A7 = 0, 



(3.10) 



(3.11) 



R.eA345 — 



1 



;Cisi{Mi~Ml) 



(3.16) 



We note that Eq. ()3.14|) can be satisfied when ai ~ 
f3 ~ 7r/4 and M2 and /i are both large compared to Mi: 



sin2ai _ M| -t- A/f 
sin 2/3 ~ Af I - 



(3.17) 



In this limit, ^ the 2HDM will correspond to the Higgs 
sector of the CP-conserving MSSM. In fact, since the 
two potentials in this limit will be the same, the other 
trilinear couplings in the two models will also be identical 
in this limit. 

We note that the MSSM couplings dSS]), hke those 
of the 2HDM, vanish for certain choices of the mixing 
angles a and (3. The exact values of parameters where 
they vanish will be modified when radiative corrections 
are taken into account. However, the modifications will 
be only quantitative in nature. 

As a special case of the limit of no CP violation dis- 
cussed above, we note that the MSSM Higgs couplings 
possess an additional important pro perty, which is usu- 
ally referred to as the decoupling [23| . This property can 
be described as follows. If Af^o ^ Mz, then a — > (3—Tr/2 
and the coupling Xhhh approaches A|^^ [28,,2§|. It fol- 
lows that in this limit, which is called decoupling limit, 
^1 1. In the present notation, this requires ai ~ /3, 
with A'/2 — Ai. (However, the correspondence between the 
2HDM and the MSSM described above, fails for tan^ 
away from unity.) We will return to the issue of decou- 
pling in Sec. IVrC] 



where g' is the U{1) hypercharge coupling. We note 
that the constraint (|3.1ip is accessible in approach (A), 
but not in (B). Using this correspondence, together with 
ai ^ a + Tr/2 and 



1 



i9'+9'')v 



gMz 



2 cos 9w 



we find 



Am ^ Xhhh, A112 ^ A/,. 



hH- 



(3.13) 



For the relations (|3.10p to be satisfied, one needs the 
non-trivial relations 

Ai=A2 = -ReA345. (3.14) 
When (PTT|) holds, and 

^2 = as = 0, (3.15) 
then (more general relations can be found in Appendix [A]) 
1 



Ai 



clMi + slMi - sl^l% 



2. Numerical study 

For fixed values of the Higgs boson masses, only a small 
domain in the paramet er sp ace of a = (ai, Q!2, 03) is com- 



(3.12) patible with positivity 2C j_ 3C |, and the perturbative uni- 



tarity in the Higgs sector 



3ll 3^.l33t. Furthermore, there 
are experimental constraints from B physics: B ~f Xs"f 



M, M\3^, B ^ TDr [37|, |38l, and B - B oscillations 
35I . l39j . These do not depend on details of the neutral 
Higgs sector, but since these processes can get contri- 
butions from exchange, they constrain the allowed 
values of tan/3 and Mfj±. Furthermore, there are exper- 
imental constraints that do depend on the neutral Higgs 
sector. These are Ri, [40l. [4H. non-observation of a light 
neutral Higgs boson at LEP I40l. and most importantly 
the constraint arising [20I |42|. |43| from Ap In con- 
trast to the case of the MSSM, {g — 2) does not play any 



1 With c2 



e, one finds to lowest order in e: c? 



e(2^j2/M| - 1) and 



2e2(M| 
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important role here. The experimental constraints on the 
Higgs sector of the 2HDM are discussed in ■ Here we 
shall follow the same approach, defining a function 



= Xgcncral + J2 ^*^(")' 



(3.18) 



where the first term, which is independent of a, is due to 
various B-physics constraints, B — s- Xg^, TVr and 

B-B oscillations 



2 _ 2 

Xgciieral 



xl 



Xb-b^ 



(3.19) 



and the second term is a sum over contributions due to 
observables Oi that depend on a. These are the non- 
observation of a light neutral Higgs boson at LEP2, the 
bb decay rate, and Ap: 



2/ N (Ci,2HDM(Q:) — C'i.rof)^ 



(3.20) 



We adopt the same definitions and experimental data as 
were used in [l7| , allowing parameters for which ^ 
5.99, corresponding to a 95% C. L. 

Some of the main features of these constraints are 
worth stressing. First of all, for small values of the 
allowed range of tan f3 is restricted by unitarity to values 
below 5-8 17[ • Secondly, for large values of Mh± , the 
Ap constraint requires M2 ^ M3 ~ Mb± , especially for 
large values tan /? . 



ki=Oj 



^1 



TREE LEVEL 



4 - 



(120,400,400,0) 
(120,300,300,0) 



tan|3=1 

a,=±jt/4, a =it/12 



(A) 




a/ (71/2) ^ 



FIG. 1: The trilinear coupling ratio ^1 defined in (|3.8|l at the 
tree level plotted as a function of ai for both the approaches 
(A) and (B), for two different values of 02 (— 7r/4, dashed, and 
+tt/4, solid), both with as = 7r/12. The masses of the Higgs 
bosons are Mi = 120 GeV, with M2 = M^± = 300 GeV 
(red, thin lines), and M2 — Mu± — 400 GeV (green, heavy 
lines). Furthermore, p, = and tan/3 = 1. The numbers in 
insets give values of {M\, M2, M^± , p) for approach (A) and 
for those of (Mi , M2, M3, Mu± , p) for approach (B). 



Within an allowed domain in the parameter space, the 
trilinear couplings can have a very strong dependence on 
the neutral Higgs boson mixing angles, as is illustrated 
in Fig. [U where we have plotted ^1 as a function of ai, 
keeping a2 = ±7r/4 and 03 = 7r/12 fixed. 

In this figure, the values of the trilinear coupling are 
calculated at tree level with Mi = 120 GeV, ^ = 0, 
tan /3 = 1, a2 = i''i'/4, and = 7r/12. The upper set 
of lines in this figure (larger values of the trilinear cou- 
plings) refer to the approach (A) of Sec. IIIBl where we 
take Ae = A7 = 0, and fix M2 = Mh± = 300 GeV, or 
M2 — Mb± — 400 GeV, as indicated in the caption. 
The lower set of lines refer to the approach (B), where 
we fix {M2,M3,Mh±) equal to (300,400,300) GeV, 
or (400,450,400) GeV, with the additional constraints 
ImAs = 0, ReAe = ReAy = 0. 

Prom this figure we note a marked qualitative differ- 
ence between the values of the trilinear coupling calcu- 
lated via the two approaches of specifying the input pa- 
rameters. In general, they both lead to solutions for only 
limited ranges of ai, keeping M3 fixed [approach (B)], 
the variation of the coupling with ai is modest, whereas 
when Ag = A7 = and M3 is calculated [approach (A)], 
the variation with ai is quite strong (upper part of the 
figure) . 

The different behavior of the trilinear coupling in the 
two approaches can be understood as follows. The trilin- 
ear couplings are linear functions of the quartic couplings 
Ai in the Higgs potential (|2.2p . In both the approaches, 
(A) and (B), the couplings Ai are rather smooth func- 
tions of ai. They are also comparable in magnitude, 
with the exception that in approach (B), where we fix 
M3, we also take as input ImAs = 0. Thus, the role 
played by Im A5 in approach (A) is played by Im Ae and 
IniAy in approach (B). This switch of role from ImAs to 
Im Ag and Im A7 has as a consequence, for the case stud- 
ied in Fig. [U that 0333 (which is quite large) switches 
sign. Since the coefficient in (|3.3p that relates 0333 to 
Am is essentially Rfr^, with R13 = sinQ;2 = il/V^, this 
sign change of the large coupling 0333 is the dominant 
effect which causes the difference between the parameter 
choices in the approaches (A) and (B) (for the set of pa- 
rameters considered). In fact, truncating Eq. (|3.3[) with 
a333 = 0, one finds the ratio of the couplings ^1 = 3.8 — 5 
for approach (A), and 3.6—3.8 for approach (B), where we 
have taken Mi = 120 GeV, M2 = 400 GeV, ^2 = -7r/4, 
and a4 = 7r/12 for both the approaches. Thus, one may 
conclude that the two different ways of introducing CP 
violation in the two Higgs doublet model lead to very 
different values for the trilinear coupling Am. 

Furthermore, we note that Re A5 depends on ai 
through M3 only 17[. Thus, in approach (A), with 
A/3 depending on ai, the resulting ReAs will also vary, 
whereas in approach (B), ReAs is constant. This con- 
tributes to the stronger variation of ^1 with ai in ap- 
proach (A). 

In view of the strong dependence of the trilinear cou- 
plings on the mixing angles, we shall henceforth average 
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over these angles within the allowed ranges. In Figs. [2] 
and [3] we show the tree- level ratios 



(A 



ml 



\SM ■ 



(6) 



(A 



112/ 



\SM ■ 
^HHH 



(3.21) 



for a representative choice of parameter sets. These are 
averages, obtained by scanning over the a ~ (ai, 0:2, as) 
space, subjecting all model points to positivity, unitar- 
ity, and the experimental constraints discussed above, 
working within the approach (A) defined in Sec. IIIBl 
and normalizing them with respect to the SM coupling, 
as defined in Eq. (|3.8|1 . The scans are performed over 
200 X 200 X 100 points. 



Before considering the general case of arbitrary CP vi- 
olation, we first consider, in Fig. [21 a region of parameter 
space near the limit p.l7p of no CP violation: 

|a2|<ao, |a3|<ao, ao = 0.05 x 7r/2, (3.22) 

referred to as "minimal CP violation" . From the above 
discussion (see Eqs. p.l0|) - p.l3|) ). it follows that this 
corresponds to a domain of parameters close to the CP- 
conserving Higgs sector of the MSSM. We note that the 
approach (A) , where we determine from Mi , M2 and 
the rotation matrix, does not permit us to go all the way 
to the limit a2 ~ 013 = 0. 



I Tree level, |i=0 | 




FIG. 2: Trilinear Higgs coupling ratios {^1) and (^2) as defined 
of the Higgs masses Mi = 120 GeV, M2 = 300 GeV and ^ = 
MSSM-like limit as defined in Eq. (|3722)l . 



I Tree level, ^=0 | 




in Eq. (|3.2ip at the tree-level, in approach (A), for the values 
0. The ratios are plotted for the case of the CP-conserving 



In this figure, which is valid for /i = 0, the plots do not extend much beyond tan/? = 5. This is caused by the 
constraint of unitarity in the Higgs-Higgs scattering sector 21]. At the higher values of tan 0, there is also a narrowing 
of the allowed region, largely due to the constraints following from Ap and the B physics [20I. [2l|. 

The surfaces representing (^1) and (^2) in Fig. [2] are in fact quite remininscent of the corresponding averages of Ai 
and Re A345 , respectively. The details of these correspondences depend on which domains of the a space (in this case, 
ai) is populated. Typically, only scattered regions are allowed [20]. 

On the other hand, in Fig. [3] we have scanned over the full range of values of ck2 and as, compatible with all the 
constraints. 

From these figures, we see that in much of the tan/? — Mfj± parameter space, there is a considerable enhancement 
of the trilinear couplings as compared with the corresponding trilinear coupling in the SM. This is consistent with 
the special case studied in Fig. [TJ 

Incorporating CP violation in the 2HDM, we see a considerable change in the behavior of the couplings as we go 
from Fig. [5] to Fig. [3] This reflects the fact seen in Fig. [1] that the couplings have a strong dependence on the mixing 
angles, and, thus, on how CP violation is incorporated in the model. 

Furthermore, there is a rapid variation of the Higgs coupling ratios around tan/3 — 0(2). As we shall see in the 
following, this is accompanied by a strong variation as one scans across the different values of a. 

The trilinear coupling has a rather complicated behaviour across the tan/?-M^± plane, despite the smoothening 
that is implicit in the averaging over a. As is evident from Fig. [T]for ^1, the ratios ^1 and ^2 of trilinear couplings 
also vary considerably over the parameter space of a. One measure of this variation is the variance, defined as 

- yOeT^Wy- (3.23) 

In Fig. [4] we display the variances ui and a2 corresponding to the upper panels of Fig. [3l This quantity (in units 
of the SM trilinear coupling) is seen to be quite considerable, in particular for ^2 and for moderate values of tan/3. 
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I Tree level, |i=0 | | Tree level, |i=0 | 




FIG. 3: Tree-level trilinear coupling ratios (^i) and (^2) as defined in Eq. (|3.2ip . in approach (A), for the values of the Higgs 
masses Mi — 120 GeV and M2 = 300 GeV. For the top panel = 0, whereas for the bottom panel fi = 200 GeV. This plot is 
for the general case including CP violation, and is compatible with general constraints on the model. 



I Tree level, |i=0 | | Tree level, |i=0 | 




FIG. 4: Variance of trilinear coupling ratios ^1 and ^2 at the tree level. The parameters here correspond to the upper panels 
(/i = 0) of Fig. El 



For approach (B), when we keep all three neutral Higgs boson masses fixed, the behaviour is smoother, as seen 
in Fig. O For the parameters considered (note that /i = 0), the values of tan/? extend up to about 5. We note 
that the quantity (^1) is rather similar to that shown in Figs. [2]and[3l increasing more or less monotonically with 
tan/3. The ratio (^2) differs more, exhibiting a bulge of rather strong H1H1H2 coupling for tan/3 = 2 — 3 and 
Mh± = 450 - 500 GeV. 
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I Tree level, |i=0 | | Tree level, |i=0 | 




FIG. 5: Trilinear coupling ratios (^i) and (^2) at the tree level [see Eq. (|3.21|l ]. for approach (B), with Afi — 120 GeV, 
M2 = 300 GeV, Ms = 400 GeV and ii = 0. 



B. Charged— Charged— Neutral Higgs Couplings 

Although in this paper we are mainly concerned with the trilinear couplings of the neutral Higgs bosons, here we 
briefly discuss the trilinear couplings of the charged Higgs boson of the 2HDM [2^ . These are given by 

Ai-i — — ^ Ri„ib„i. (3.24) 

m=l,2,3 

We note that the index m (like the indices n and o in Eq. (j3.3p ) refers to the weak interaction eigenstates (rji, 7]2 and 
773) of Eq. (|2.4|) . The coefficients in (|3.24p can be written as (in units of —iv): 

hi = cos/3{sin^ [3 (Ai — A4 — Re A5) + cos^ /3 A3 
+ cos P sin/3[(tan2/3- 2)ReA6 +RCA7]}, 

62 — sin/3{cos^ [3 (A2 — A4 — Re A5) + sin^ [3 A3 
+ cos P sin /3 [Re As + (cot^ /3 - 2) Re A7]}, 

63 — cos/3 sin/3 ImAs — sit? film \q — cos^ /3Im A7. (3.25) 

In Fig. [5] we display the coupling of Hi to the charged H^H~ pair, for the same choice of parameters as in Fig. [T] 
We note that in approach (B), with tan/3 = 1 and ImAs — 0, this coupling has two interesting features: (1) It does 
not depend on the sign of 02; and (2) it passes through zero for ai = — 7r/4. The independence from the sign of a2 
can be understood in the following manner. For i = 1, the first two terms in Eq. p.24p (proportional to -Rn and R12) 
contribute terms which are odd in a2, arising from Ai, A2 and A3. These odd parts are proportional to (M| — M|), 
and cancel against odd parts coming from the third term in (|3.24p . which are proportional to ImAg + ImA7, which in 
turn is proportional to 7W13 + 7U23. 

Furthermore, the vanishing of the coupling at ai — — 7r/4 can be understood as follows. In approach (B), with 
tan/3 = 1 and ai = — 7r/4, we have (see Eqs. (3.1) and (3.2) of Ref. [iTj ) 

Ai(a2,a3) = \[clMl + (C3 - S2S3)'M| 

+ (s3 + S2C3)'A/|-m2]^ 

A2(a2, 03) = \[clMl + (C3 + S2S3)^Ml 

+ (s3-S2C3)2a/2-/i2]. (3.26) 

Thus, for these particular values of tan/3 and ai, A2(— a2iQ^3) = ^i{<^2t cts), the first two terms in Eq. (I3.24p cancel, 
and we are left with 



Ai+_ = Risbs = is2(ImA5 - ImAe - ImA7) = 0, 



(3.27) 
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where in the last step we have used Im A5 = 0, and the relations (see Eq. (|2.19p ) 



/2 

Im A5 + Im Ag + Im Ay = I-^ 13 + -^23] . 



(3.28) 



valid for tan /3 = 1 , and 



(3.29) 



valid for ai — — 7r/4. 



Ui=o 



Charged Higgs Coupling, TREE LEVEL 




FIG. 6: The trilinear coupling Ai_|._ at the tree level as a function of ai, for two different values of 02 (— 7r/4, dashed, and 
+7r/4, solid), both with 03 — n/12. The masses and parameters are the same as in Fig. [l] 

In Fig. [71 we show the coupling of the lightest neutral Higgs bosons to the charged Higgs bosons, averaged over a, 
in approach (A). Corresponding to what we observed in Fig. [6l these averages over a also vanish along certain lines 
in the parameter space. 





450 
% 400 
'^eiy 350 



FIG. 7: Tree-level trilinear couplings Ai+- and A2+- (in units of —iv) of the two lightest neutral Higgs bosons. Hi and H2, to 
a charged pair for Mi — 120 GeV, M2 = 300 GeV and fj, — 0. The couplings are averaged over a, and have been plotted for 
approach (A). 



IV. ONE-LOOP CORRECTIONS TO THE loop radiative corrections to the tree-level potential 

HIGGS BOSON MASSES When this is done, the relations among masses and the 

A. One-loop effective potential 



For a realistic computation of the Higgs boson prop- 
erties in the 2HDM, we must take into account the one- 
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couplings of the tree-level potential undergo a change. 
In the MSSM, this effect, which is dominated by the 
top (and stop) quark contributions, is known to be very 
important [lij (see also Ref. ^). In the 2HDM, there 
being no one- loop contribution from squarks, the main 
contributions at one loop are generated by the Higgs 
bosons. We shall here study the one-loop radiative cor- 
rections to the neutral Higgs sector using the method 
of one-loop effective potential [13] ■ The one-loop correc- 
tions to the tree-level potential (|2.2p are given by 



647r2 

E 

form ions 



E 

bosons 



, M2 3 



f2 



log 



M 



(4.1) 



where the sums run over all bosonic and fermionic degrees 
of freedom, respectively, and represents the scale at 
which the couplings are evaluated. As indicated above, 
the dominant contributions to (|4.ip come from the Higgs 
bosons (neutral and charged), and the top quarks, all 
of which we will calculate in this section. Taking into 
account the spin, charge and color degrees of freedom, 
we obtain 



Ay 



^^^cutral Higgs 

1 



charged Higgs 



A1/t 



top 



3 

647r2 E 

1=1 



/r4 



, A/2 3 
log 



log 



167r2 



log 



(4.2) 



adjusting the parameters of the potential such that the 
minimum is stationary in terms of the same vacuum ex- 
pectation values vi and V2 as for the tree-level potential 
(i.e., tan/? remains unchanged). In the following we shall 
calculate the different contributions to the neutral Higgs 
boson masses arising from ()4.2p using the above proce- 
dure. 



B. One-loop corrections to Higgs masses 

The mass squared matrix for the neutral Higgs bosons 
will receive contributions from one-loop effects due to 
neutral Higgs bosons, which can be written as 



d AT^cutral Higgs 

drjidrij 



1 



3 r 



327r^ 



E 



(32 Q^^ Qjj. 



+ Mf(log^-lj^ 



(4.5) 



The derivatives of Mf with respect to the fields rji, 
which are rather involved, are discussed in detail in Ap- 
pendix [cl 

The corresponding second derivatives of the charged- 
Higgs-field-dependent part of the one-loop corrections to 
the potential, Aychargcd Higgs, in (|4.2p are given by 



^ AV^}^ai-gc(i Higgs 

driidrij 



1 



167r2 



(32 5^. Qjj. 

drjidrij 



M|f±(log-^^^* 



1 



(4.6) 



where all masses are understood to be field-dependent 
masses. Thus, after differentiation analogous to the one 
in Eq. (|2.9p for the tree-level potential, for example, we 
do not set the fields 771, 772 and 773 to zero, but instead 
evaluate 

MUm,m.Vi) - 7^7^n'7i,^2,%,i?^,G±), (4.3) 



and similarly 



Mh± (771, '72, '73) 



_d d_ 

dH+dlf- 



(4.4) 

and after having performed the differentiation, we set 
ff± = G± = 0. 

At tree level, the three minimization conditions (|2.18p 
(two real, and one imaginary) can be used to eliminate 
the bilinear parameters 7^2 and ttt.^, as well as relate 
ImAs to Im77^f2• At the one- loop level, the expressions 
for 777,2 gj^^ ^2 g^j.g niodified due to the contributions 
coming from dAV/drji and dAV/dri2- This amounts to 



Finally, there are one-loop corrections due to the top 
quark, which can be written as 



where in the last equation we have used dmt/drji = 
(rnt/v2)Si2. We neglect the contributions from lighter 
fermions. 

The full one-loop mass-squared matrix for the neutral 
Higgs bosons can, thus, be written as a sum of four terms: 



-^Lc,jj(m,?72,??3) 



^V^cutial Higgs 



C^^AVcharged Higgs 

1- — • (,4.8) 



dr]idr]j 



dr]idrij 



From the mass-squared matrix (|4.8p . the eigenvalues 
Ml (771, 772, 773) can be determined. These are the one- loop 
corrected masses for the neutral Higgs bosons. 
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In Fig. [8] we show, for approach (A), the one- loop corrected Higgs boson masses, corresponding to the trihnear 
couphngs studied in Fig. [1] (i.e., with a2 and held fixed), due to the one-loop radiative corrections discussed above. 
The superscript 'LO' here denotes lowest-order values. The Higgs boson masses are determined as eigenvalues of the 
matrix (|4.8p . evaluated for 771 = 772 = % = 0. The full 77-dependence, to be discussed in the next section, will however 
be required for evaluating the one- loop trilinear Higgs couplings. 

Approach (A) 

600 
500 
400 
300 
200 



100 

-1 -0.5 0.5 1 

a, / (71/2) 

FIG. 8: One-loop Higgs masses Mi, M2 and M3 (in units of GeV) plotted as functions of qi, with Q2 = ±7i"/4, both with 
a-i = 7r/12. The lowest-order Higgs masses are Mi^° = 120 GeV, with = M^jl = 300 GeV (thin lines), and = 

M]^ = 400 GeV (heavy lines). Here ji — Q and tan/3 = 1. Dashed lines refer to tree- level values, whereas green and blue lines 
refer to one-loop-corrected results, with the scale parameter Q — 2Mi'~' and Q = 500 GeV, respectively. 

The following points are worth noting here: (i) The tree-level value of M3, namely Ml"*-*, has a strong dependence 
on ai. This is a consequence of keeping M]"^, Af^"'-', a2 and fixed, (ii) The one-loop corrections have a significant 
dependence on the scale at which they are evaluated, as is illustrated by comparing the values of the Higgs boson 
masses at two different values of . (iii) The one- loop correction to the Higgs boson mass can be either positive or 
negative, (iv) At some of the boundaries of the allowed parameter space, these one-loop corrections become large. 

The reason for the strong variation of il/3 at the edges of the allowed parameter space can be found from an 
inspection of the different origin of such edges. These can be described as follows: 

1. In approach (A), at tree level, (M^*-*)^ is determined as a rational fraction, where the numerator is linear in 
(M^f"*-*)^ and (M|"'^)^ (multiplying rotation matrix elements and tan/3), whereas the denominator is a linear 
function in rotation matrix elements and tan/3 (For an explicit formula, see Eq. (4.16) of Ref. jl^.). For critical 
values of the rotation matrix elements, e.g., when R31 = R32 tan/3, the denominator will pass through zero. At 
such points, M^'^ diverges, and the allowed parameter space terminates. 

2. In approach (A), Af| may at some critical values of the rotation matrix elements become equal to Af|. This 
happens when either i?i3 = or i?i2tan/3 = Rn [Toj . In the former case we have a2 — (ai and being 
arbitrary), whereas in the latter we have tan/3 = cotai (a2 and aa being arbitrary). 

3. Positivity may break down. 

4. Perturbative unitarity may break down. 

5. Parameter space may be truncated by some experimental constraint. 

The first two of these cases are illustrated in Fig.[51 We actually exclude, as 'unphysical', regions where Mi < 100 GeV 
(for M^^ = 120 GeV) and where either M^'=' > 2M^'^ or M3 > 2M2. 

For approach (A), positivity is expressed in terms of simple inequalities, whereas for approach (B), with and 
Qfy non-zero, an involved numerical test is required [l^]. Perturbative unitarity in the Higgs-Higgs scattering sector 
tends to truncate high values of tan/3, unless ^ is large (comparable to Af2) [l7i]. 

Having described the radiative corrections to the Higgs masses, we next show in Fig.lHltypical one- loop multiplicative 
mass corrections, i.e., Mi/M^''^ and Mz/A/^^o for the case Af|^0 = 120 GeV, M^'^ = 300 GeV, and ^l = 0. (This will 




14 




be discussed in more detail in Sec. IVIin The renormalization scale has been chosen as = {2M\'^Y . As in Section 
mil an average over allowed points in a has been performed. 

For Ml, the loop correction is (for = (2AfJ"'-*)^) negative at low values of tan/3, whereas it is positive for 
higher values. For Af2, on the other hand, the correction is smaller, and has no strong dependence on tan/?. The 
negative correction for Mi is different from the MSSM, where the one-loop corrections, dominated by the top-quark 
(squark) contribution, is positive. Here, it is the Higgs sector which dominates, and this contribution has the opposite 
sign, because of the bosonic nature of the contributing particles in the loop. In actual practice the situation is 
more complicated, since different Higgs fields may contribute with different signs, depending on how the mass of the 
respective Higgs boson compares with the scale Q. 

In view of the fact that the one-loop correction to Mi can be significant, and negative, we impose the constraint of 
LEP2 non-discovery of a Higgs boson on the one-loop corrected lightest Higgs boson mass. Higgs boson masses below 
the "magic" 114.4 GeV value are of course allowed, provided the coupling to the Z and bb (or r pairs) are suppressed. 

Treating the loop corrections as a perturbation, we have here imposed the theoretical and experimental constraints 
on the tree-level parameters. We note that at low values of tan/3, the mass of the lightest Higgs boson. Mi, may be 
considerably reduced. When we proceed to study the trilinear couplings, we shall therefore impose the LEP2 non- 
discovery constraint on the loop-corrected mass, Mi. For a higher scale parameter, though, the mass may increase 
also at low values of tan (3. 



V. ONE-LOOP CORRECTIONS TO THE TRILINEAR HIGGS COUPLINGS 



The one-loop corrected trilinear Higgs couplings are given by the third order derivatives of the one-loop effective 
potential with respect to different fields. We shall separately consider the contributions to the potential from the 
neutral Higgs fields, the charged Higgs fields and the top quark fields. We start the calculation by evaluating the 
derivatives with respect to the weak fields, which can then be converted to derivatives with respect to the fields 
corresponding to the mass eigenstates with the help of Eq. p.2p . 

The starting point of the calculation of the one-loop corrections to the trilinear Higgs couplings is the expression 
(|4.2p . From this we obtain the one-loop radiative corrections to the trilinear Higgs couplings due to neutral Higgs 
bosons as 



neutral Higgs 



1 



dr],dr]jdr]k 



32n' 



E 



1 dM'^dMjdMj 
Ml drji drjj drjk 



log 



Mj 



d^Mf dMj ^ d^Mf dMj 



\driidrij drjk 
d^Mj dM? , 

IVif 



dilkdr]i dr]j 



drjjdr]k d-q, 
M? 

Q 



1 



dvidmdiik 



(5.1) 



The derivatives of squared masses with respect to the weak fields iji in (jS.ip are evaluated and discussed in Appendix [Cl 
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To find one-loop charged Higgs boson correction to the trilinear couplings, we start with field- dependent squared 
masses, (771, 7721 %)i for the charged Higgs boson. Noting that the potential is a function of masses, i.e., 

^ Charged Higgs — ^y(,Vij''l2,V3), we obtain 



harecd Higgs 



d-q^drjjdrik 
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log 



Ml. 



dHll^ dMl^ 
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(5.2) 



For the one-loop corrections arising from the top-quark, we note that in Model II (which we consider) the up-type 
quarks couple only to <i>2 (and not to $1). Thus, the one-loop contribution due to top quark affects only the fields 
772. We can then write 



d^AVt 



top 



12 mf 



drjtdrjjdijk tt^ 



1 '^t c r r 
Iog-rj0j2'5j20fc2- 



(5.3) 



VI. QUANTITATIVE RESULTS FOR TRILINEAR HIGGS COUPLINGS 



We have seen that the lowest order trilinear couplings are highly dependent on the details of the two Higgs doublet 
model. This sensitivity on the details of the model is actually enhanced when we take into account the one-loop 
corrections to the Higgs potential. The loop corrections to the trilinear couplings can vary strongly with the parameters 
of the model, and also have a considerable dependence on the scale used in Eqs. (|5.ip - (|5.3p . 



A. Different contributions 



For fixed values of 02 and as, as in Fig. [U we show in Fig. [TO]the ratio of the one-loop corrected trilinear Higgs 
coupling Am and the corresponding tree-level coupling. This ratio is highly dependent on ai (as is the coupling 
itself). 
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FIG. 10: Contributions to the trilinear coupling ratio plotted as functions of ai, for 02 ~ ±7r/4, both with as — ■k/12. The 
Higgs mass and parameters are as in Fig. [1] The scale parameter is Q = 2Mi"°. The full correction is shown by the upper, 
heavy, black lines, whereas partial contributions from neutral and charged Higgs bosons, as well as from top quarks, are shown 
by the thin, red, green and blue lines. 
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As compared with the one-loop corrections to the masses, the corrections to the trihnear couphngs are relatively 
much larger. The following two points are worth noting: (i) Most of the one-loop contribution comes from the 
neutral-Higgs term, Eq. (|5.ip (red lines labeled "ffi, H2, ffa"). (ii) The correction can be very large at the boundary 
of the physically allowed region (in this case, in ai). These regions of very large corrections, which are excluded from 
averages over a, will be discussed further in the following. 

In approach (A), when we scan over the values of a, there are regions where M3 gets very close to M2. From 
Ref. [19], we see that (at tree level) this happens when (i) 02 — > or (ii) tan/3 cotai, or (iii) a2 ±7r/2. Case 
(ii) occurs for tan/3 = 1 and ai = tt/A. In this region, the trilinear coupling tends to get very large. This is clearly a 
region of parameter space where the present approach breaks down. In the notation of Appendix [Cl [see Eq. (jC12p ]. 
we see that two roots merge (M3 M^) when the complex phases of Si and S2 differ by 27r/3. Near such points, the 
derivatives of Si and S2 with respect to a and &, particularly the second and third-order ones, tend to get very large. 
Actually, such points have been removed from the following figures, and we impose a cut when scanning over a in 
Sec. lvfDl 

uim{{Ml° ~ M2^°), (A/3 - A/2)} > 5 GeV. (6.1) 
B. Scale dependence 

Just as in the case of one-loop corrections to the masses, we note a strong dependence of the trilinear couplings on 
the scale parameter that enters the effective potential. This scale plays a role somewhat analogous to the mass scale 
of superpartners in the effective-potential approach to the MSSM. 




FIG. 11: Trilinear coupling ratio [see Eq. (I3.8|l ] as a function of ai. Colored lines refer to tree-level results (same color 
coding and parameters as in Fig. [T]), black lines refer to loop-corrected results (thin: Q = 2MJ"°, heavy: Q — 500 GeV). 

In Fig. [TT]we study the one-loop corrected trihnear couplings, for both the approaches (A) and (B), and for two 
sets of Higgs masses within each approach (the same as in Fig. [T|) . We show how the one-loop corrections depend on 
the scale Q. For this purpose we have considered two timelike values of Q, namely, Q — 2A-/|"° and Q — 500 GeV. As 
in the case of Fig. [TOl we see that in approach (A), the couplings tend to get large near the boundary of the allowed 
parameter space. (Near ai = 7r/4, we notice the irregular behavior discussed in Sec. IVI AI when A/3 Ah-) 

The dependence on the scale is rather involved, but there is a tendency that a lower scale gives a larger correction. 
This is rather clear in the upper part of the figure, which is a plot for approach (A). For approach (B), however, we 
have not labeled the individual curves, the main point being to show that there is a considerable uncertainty related 
to the choice of scale in both approaches. 

C. Decoupling 

In the context of the trilinear Higgs coupling, decoupling would imply that the self-coupling of the light Higgs boson 
could be expressed in terms of its mass, as in the Standard Model. This was indeed found to hold for the MSSM, if 
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FIG. 12: Trilinear coupling ratio (^i) as defined by Eq. (|6.2|l . Diflerent quantities in approach (A) with Mi'^ = 120 GeV, 
= A* = 500 GeV, and Q = 2MJ'°. The ratio (^i) ~ 1 demonstrates approximate decoupling for M2 = p ~ Mjj±. 



In order to discuss decoupling in the 2HDM, we first need to define, as a reference, a Standard-Model-like limit by 
imposing the tree- level constraint (|3.1ip . together with Eq. (I3.22p . The SM-like nature of this limit may be extended 
to the loop level, if in the loop corrections (|4.2p we include only the contributions of the lightest neutral Higgs boson 
Hi and the top quark. Within this framework, we may then study decoupling, by defining the ratio, analogous to 



a full \ 

(Ci) - (6.2) 



where 

• A™ includes only loop corrections due to tt and Hi, as would be the case in the SM. 

• X^^il includes all one-loop corrections to the Higgs coupling in the two Higgs doublet model, i.e., also those due 
to H2, H3 and 

• The averaging over a2 and is constrained by Eq. (|3.22p with uq — 0.025x7r/2, whereas ai is left unconstrained. 

It should be noted that here we do not insist on a correspondence between the 2HDM and the MSSM, since this 
would require also ai ~ p. 

If the ratio (|6.2p were to come out as unity, it would be a signal for decoupling: the heavy Higgs bosons do not 
affect the coupling of the lightest one with itself. In Fig. [121 we display this ratio, for rather high values of M^"*-* and 
/x: M'l'^ =11 = 500 GeV, subjecting all parameters to the general constraints discussed in Sec. lIII Al For A/^± ~ Af2, 
the ratio (|6.2p is indeed rather close to unity, hut for larger values of Mfj± , it becomes significantly larger. In the 
region where (^1) ^ 1.20, we have also checked that M3 — M2 < 0(40) GeV. Thus, in this region of the parameter 
space, decoupling holds at the loop level to about 10-20%. For larger values of Af^±, however, where the model is 
still consistent, decoupling is strongly violated. 



The result displayed in Fig. [12] is quite stable under 
a change of ao from 0.025 x 7r/2 to 0.05 x n/2. How- 
ever, the result depends on the choice of scale adopted 
for evaluating the one-loop potential. For example, for 
Q = 500 GeV, (^1) differs significantly from unity. This 
is perhaps contrary to expectations, since log(Af|''-'/(3) 
then vanishes. However, we recall that (i) the one-loop 



potential also contains non-logarithmic terms, and (ii) 
the one-loop potential modifies the minimization condi- 
tions, i.e., its contribution shifts the values of the soft 
parameters m^^ and 77132 ^f the tree level potential. 

In the 2HDM, in contrast to the MSSM, the decoupHng 
cannot be exact. The reason is obvious: the heavy fields 
that we would like to neglect, are all bosonic, and thus 
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I Loop level, |i=0 | | Loop level, |i=0 | 




ratios (^i) and (^2) [see Eq. (|3.2ip ]. Middle: one-loop multiplicative mass correction. Bottom: Variance of trilinear coupling 
ratios ^1 and ^2- 



all contribute with the same sign to the potential. Fur- 
thermore, even if we choose the scale Q such that the log- 
arithms vanish, there would still be the non-logarithmic 
remainders. 

The question of decoupling can be illustrated by the 
behavior of the p parameter, which gets contributions 
from all pairs of Higgs bosons. These contributions van- 
ish for equal masses, but are large when the masses are 
very different. 

In the limit of no CP violation, custodial symmetry [i^ 



would imply M^± = A/j, where Mj is either the mass 
of the CP odd {A) or a CP even {H) Higgs particle 
In the case that M^± = A, orthogonality of the even- 
sector mixing matrix would protect the Ap parameter 
from large corrections. In the present case of CP viola- 
tion, a higher degeneracy is required, M2 = A/3 = Mfj± , 
in order to avoid large contributions to Ap. Since this is 
in general not satisfied, the experimental constraints on 
Ap severely constrain the allowed parameter space [2l| . 
and there is no decoupling. 



19 



I Loop level, n=200 GeV | 




I Loop level, t^=200 GeV | 




FIG. 14: Same as Fig. [TS] 



D. The general case 

After this exploratory discussion, we are now ready to 
scan and average over the full range of a, compatible with 
theoretical and experimental constraints on the model, as 
was done at tree level in Sec. IIII Al Figs. [T3] and [H] are 
plotted for approach (A) with A/^o = 120 GeV, M^'^ = 
300 GeV, = (Fig. [TSl), a* = 200 GeV (Fig. [HI) and 
with scale Q = 2M|"° in both cases. 

In these figures, we show the ratios (^i) and (^2) of 
trilinear couplings in the top panels. The middle panels 
are devoted to the loop-induced corrections to the masses 



I Loop level, n=200 GeV | 




I Loop level, |x=200 GeV | 




but with fi = 200 GeV. 

I 

Ml and M2, whereas the bottom panels show the vari- 
ance associated with the averaging of ^1 and ^2 over the 
a parameter space. 

An important difference between the two cases pre- 
sented in these figures is the range in tan f3. As already 
discussed, for a low value of /x, the high values of tan (3 
would lead to a model that would violate unitarity, and 
are, therefore, excluded [20j . 

Comparing the top panels of Figs. fT3l and [Ml with the 
corresponding tree-level resuls shown in Fig. [31 we note 
the following: (i) The overall shape in the tan/3-Mj:^± 
plane is similar, (ii) The loop-corrected couplings are 
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somewhat larger, after averaging over a, in analogy with 
the results shown in Fig. [11] for a2 and fixed. 

The center panel of Fig. [13] is rather similar to Fig. [9] 
the only difference being that here the LEP2 constraint 
is imposed at the loop level. Focusing first on the loop 
corrections to Mi, and comparing the mass correction 
at the two different values of fi, we note that the case 
jj, = gives a larger reduction at low values of tan (3 and 
a larger increase at high values of tan {3. For M2, on the 
other hand, the loop correction has only a rather weak 
dependence on tan/3 (and Mh±)- For a higher value of 
the scale Q, however, the loop correction to Mi may be 
positive also at low values of tan f3. 

The trilinear couplings represented in these figures, are 



of course not physical, they are instead "typical values" , 
obtained by an averaging over the angles of the mixing 
matrix R that diagonalizes the mass-squared matrix. As 
discussed above, the variation over the a parameter space 
can be considerable, as illustrated by the variance shown 
in the bottom panels. 

A more detailed comparison reveals that the loop cor- 
rections are most important for low values of low-tan /3. 
This is particularly so in the case of (^2), where the vari- 
ance is also quite large. This is in part caused by the 
ill-defined limit M^ M2 (02 0) discussed above. 
[Actually, at high values of M[j±, the (^2) plots are trun- 
cated at the top.] 



I Loop level, |i=500 GeV| | Loop level, 1^=500 GeV| 




FIG. 15: Similar to the upper panels of Figs. [13] and [14] but for /i = 500 GeV. Upper panels: Mf^ = 300 GeV; lower panels: 
M^"'-' — 500 GeV. Note that the scales are different from those in the earlier figures. 



The correlation between these trilinear couplings and 
the quartic couplings Xi of the tree-level potential or their 
averages, (Xi), is not very simple, but a couple of observa- 
tions can be made. In Fig. 1131 for n = the enhancement 
at low values of tan P is correlated with a corresponding 
enhancement of (A2), whereas (Ai) is large at high values 
of tan (3. 

For the Hi Hi Hi coupling, represented by (^1), we note 



that for /X = it is rather large for all values of tan P and 
Mh± where the model is consistent, whereas for /i = 
200 GeV it falls to zero for large values of tan /3. 

Finally, in Fig.[T3]we show the behavior for somewhat 
higher value of the soft mass parameter, /i — 500 GeV. 
We display the averages (^1) and (^2) for two values of 
M^"*-*, 300 GeV (as was adopted in the earlier plots) 
and 500 GeV. Compared to the corresponding panels of 
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Figs. [Island [Ml the present ones (for = 500 GeV) are 
more smooth, and also somewhat different in magnitude. 



VII. SUMMARY AND CONCLUSIONS 

The two Higgs doublet model is a viable extension of 
the Standard Model Higgs sector. Apart from having 
a richer structure of Higgs bosons and Higgs couplings, 
it can accommodate CP violation beyond that of the 
Kobayashi-Maskawa mechanism in the SM. In this pa- 
per we have studied in detail the trilinear couplings of the 
lightest Higgs boson of this model, including the CP vio- 
lating effects. Our main objective has been to study the 
implications of CP violation for the trilinear couplings 
of the general two Higgs doublet model. We have dis- 
cussed two different parametrizations of the general two 
Higgs doublet model, and imposed the constraints on the 
parameters of the model following from theoretical con- 
siderations and from different experiments. Tree-level 
unitarity, B physics and Ap are the most constraining 
ones for the parameters of the two Higgs doublet moddel 
with CP violation. 

Within the allowed domain of the parameter space of 
the model, the trilinear couplings can have a very strong 
dependence on the neutral Higgs boson mixing angles. 
Because of this strong dependence of the trilinear cou- 
plings on the mixing angles, we have presented averages 
over the allowed ranges for these angles. Together with 
the variances, these show the range of trilinear couplings 
that can be expected in the two Higgs doublet model 
with CP violation. 

There is a region of parameter space of the 2HDM, 
which we call ^^minimal CP violation,^^ which corresponds 
to the domain of parameters close to the Higgs sector 
of the CP-conserving minimal supersymmetric standard 
model. In much of this parameter space there is a consid- 
erable enhancement of the trilinear couplings of the two 
Higgs doublet model as compared with the corresponding 
trilinear coupling of the SM. 

We have computed the radiative corrections to the 
masses of the Higgs bosons and the trilinear couplings in 
the one-loop effective potential approximation. The one- 
loop corrected mass of the lightest neutral Higgs boson 
can be modified by as much as 20%, whereas the heavier 
ones are less affected by the radiative corrections. The 
one-loop corrected trilinear couplings are typically en- 
hanced by the loop effects, but so is the sensitivity to 
the neutral-sector mixing angles. A similar enhancement 
due to CP violating effects has also been found for the 
MSSM [ISIIH. 

To a good approximation, decoupling holds in a limited 
range of parameters, where the model is close to the CP- 
conserving MSSM. In general, however, the heavier states 
also have a considerable impact on the properties of the 
lightest Higgs boson. In this sense there is no decoupling 
in the two Higgs doublet model with CP violation. 
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APPENDIX A: QUARTIC COUPLINGS 
EXPRESSED IN TERMS OF MASSES 

In the 2HDM, with ImAs ^ 0, the two masses Mi 
and M2, together with a and tan/3, determine M3. The 
masses will also be related to the quartic couplings of 
the Higgs potential, via linear relations that can be un- 
ambiguously inverted. It is convenient to distinguish two 
cases, whether or not Ag and A7 are non-zero. 



1. Approach (A), Ae = A7 = 

This case is referred to as approach (A) in the text. 
Providing also M^± and in addition to Mi, M2 and 
a, all the A's are determined as follows (TtI. [20j: 

+ (C1S2C3 - siS:ifMl - slfi^], (Al) 
A2 = -^[slcjM^ + (cica - siS2S3)^M| 

+ (ci.53 + sis^c^ fMl - cIp\ (A2) 
A3 = —^{ciSi[clMl + {slsl - cl)Ml 

+ (S2C3 - sl)Ml] + S2C3S3(c? - s\){Ml - Mi)} 

+ ^[2Af^± -m'], (A3) 

A4 = ^[slMl + clslMi + clclMl + - 2Af^±], 

(A4) 

Re As = 4 [-^'Mf - clslMi - cjcjMl + p^] , (A5) 

ImAs = ^{ci3[ciC2S2Mf - C2S3(ciS2S3 + SiC3)M| 

+ C2C3(siS3 - CiS2C3)M^] + Sf3[siC2S2M^ (A6) 
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+ C2S3(ciC3-SlS2S3)Af"|-C2C3(ciS3 + SiS2C3)M|]}, 

where — cos/3, sp = sin/3. 

While M| is given in terms of , M|, i? and tan/3, 
it is more transparent not to substitute for M| in the ex- 
pressions (|All) - (jA6p . These equations are the analogues 
of those of ^53] for the CP-conserving 2HDM. 



2. Approach (B), Ae 7^ 0, A7 7^ 

In this case, referred to as approach (B) in the text, 
we take Mi, M2, M2, a, Mff± and together with 
ImAs, ReAe and ReAy, as the input. In order to keep 
the notation compact, it is convenient to introduce the 
following abbreviations: 

Re A345 = A3 + A4 + Re As , ( A7) 

Re A567 = Re A5 + cot /3 Re Ag + tan /3 Re Ay . ( A8) 

Then, the A's can be determined from the following rela- 
tions: 



1. The trihnear coupling Am 



Ai - 



1 



7^ — (3c,ReA6-s, 



A2 — —T 



MI2- cjfi^ 



^ReAr) 
(A9) 

-^(clReAe-S^ReAr) 
(AlO) 



Re A; 



1 



•345 

A4 
Re As 

Im Ag 
Im A7 



3 
2 

V 



^(c^RcAe + s^R-eAy) 
Jfi^ - Mfj±]-ReX. 



567, 



(c^Re Ag -I- s^Re A7), 



1 
1 

'2^ 



2M?, 



2Xi3 



s^Im As 



c^ImAs 



(All) 
(A12) 

(A13) 
(A14) 

(A15) 



Invoking Eq. (|2.14|1 . the A^f^ can be expressed in terms 
of R and the masses Mi, M2 and M^^. 



We start by writing down Eq. p.3p explicitly for 
{i,j,fc} = {1,1,1}: 



Am— RlmRlnRlo O-mno 

m<n<o— 1,2,3 

~ 3!{i?ii[i?ii am -1- R12 0112 -f- i?i3 0113] 



-I- i?ii[i?i2ai22 + R12R13 ai23 + Rf^ 0133] 

+ Rl2[Rl2 0,222 + R12R13 ^223 + ^13 ^233] 

+ i??3«333}, (Bl) 

where the factor 3! is due to the fact that we here have 
couplings to three identical Hi fields. Substituting for the 
elements of the rotation matrix, Eq. (I2.13P , one obtains 

Am = 3!{(ciC2)^[ciC2 am -f- S1C2 aii2 -I- S2 0113] 
+ ciC2[(siC2)^ ai22 + S1C2S2 0123 + si 0133] 

+ SiC2[(siC2)^ 0222 + S1C2S2 0223 + 0233] 

+ S2a333}. (B2) 

Let us now consider the case of Ag = A7 = 0, and 
substitute for the coefficients Umno from Eq. (|3.4p . In 
this case. 

Am = 3{ciC2C^(ciC2 -I- sls'p)Xi + SiC2Sp{slcl + slcp)X2 

- 2c2sl{cici3 + siS/3)Re A34S 

— S2[(c2 — S2)c/3S/3 -|- 2ciSiC2]Im As}. (B3) 

In the limit of no CP violation, with Hi being odd, we 
have C2 = [see Eq. (|2.15p ]. and the trilinear couphng 
simplifies further to 



Am = 3s2C;3S,3lmA5 



(B4) 



with S2 — ±1. Similarly, in the limit of no CP violation, 
with H3 being odd, we have S2 — S3 = 0, and obtain the 
simple expression given in Eq. p.6ap . whereas the limit 
of no CP violation with H2 being odd, leads to no further 
simplification of the expression (|B3p . 



2. The trihnear couphng A112 



We shall again first write Eq. (|3.3p explicitly for 
{i,j,fc} = {1,1,2}: 



APPENDIX B: TRILINEAR COUPLING FOR 
SOME SPECIAL CASES 

In order to get a feeling for the trilinear Higgs couplings 
of the 2IIDM, we shall here explicitly write down the two 
trilinear couplings Am and A112 for some special cases. 



A112 — ^ ^ Ri' mRj' nRk' o O^mno 

m<7i<o— 1,2,3 

~ 2 {RlmRlnR2o + RlmR2nRlo 

m<7i<o— 1,2,3 

~t~ R2inRlnRlo\ O,rnno 

= 2{3i?iii?2i am -f Rii{RiiR22 + 2i?i2i?2i)aii2 
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-Rii(-Rii-R23 4 

-Rl2(2i?lli?22 
[i?ll(i?12i?23 - 
i?13(2i?iii?23 
i?12[i?12^23 + 
^13(2i?12-R23 



2i?i3i?2l)oil3 
f i?12-R2l)ai22 

- R13R22) + ^12-Rl3-R2l]ai23 
f i?13i?2l)ai33 + 3i?j2-R22 ^222 
2i?i3i?22]a223 

f i?l3-R22)a233 + 3-Ri3i?23 0333}- 



(B5) 



If we now substitute for the elements of the rotation ma- 
trix, for the coefficients amno, and set Ag = A7 = 0, we 
obtain the result quoted in Eq. (|3.6bjl . 



APPENDIX C: DERIVATIVES OF Mf 



The higher derivatives can likewise be written as 



dij^drijdrik 



EE 



daadafi drji drjj 



daa drjidrij 
(C7) 



E EE 

Q ^ /3 '- 7 



daa dap da^ 



daadapdaj drji drjj drjk 



d'^Mg / d'^tta dap d'^Ua dap 



daadap ydrjidrij drjk drjj drjk drji 



d'^aa dap 
drjkdrji drj. 



dM^ d^a^ 
daa drj^drjjdrjk 



(C8) 



The loop corrections to the trilinear couplings depend 
on derivatives of the squared masses with respect to 
the weak fields rji. In order to obtain these derivatives, 
a considerable amount of book-keeping is required. The 
complication arises from the fact that we allow for CP 
non-conservation and Mf will thus be determined by the 
roots of a cubic equation rather than a quadratic one. 
We start by considering the M| as functions of the coef- 
ficients of the corresponding cubic eigenvalue equation: 



dMf 



where M? 



drji 

Xf is a solution of 



d 

— Mfia,b, c) 
drj^ 



\^ + a\^ + h\ + c = 0, 
with field-dependent coefficients: 



(CI) 



(C2) 



a ^ a{rii,Tj2,rj3), 6 = 6(771,772,773), c = 0(771, 772, 7^3). 

(C3) 

These coefficients of the cubic eigenvalue equation are 
obtained from the derivatives (|2.9p . where, in contrast to 
Eq. (|2.19p . we do not set the fields to zero. In terms of 
the 3x3 mass squared matrix A^^, these coefficients of 
the cubic equation are given by 

a = -TtM^, 

6=i{(TrM2)2_Tr[(A^2)2]|^ 
c=-detA^^ (C4) 



The derivatives (jCip can thus be split up into simpler 
entities: 

— M2(a 6 c) - ^Ml^ + dMldb_ ^ dMf dc 
drji ' ' da drji db drji dc drji 

= Ef^?^' 

^ daa drji 

a 

where we have introduced the collective notation 

Oq. = {a, 6, c}. (C6) 



1. Cubic equation 

In order to obtain the derivatives of with respect to 
the aa that enter Eqs. ((C5|l . ((C7| and ((C8)) . we start by 
solving the cubic equation (jC2p in terms of the following 
notation. Let 



(7=— 6 a, r — —iab — 6c) a 

^39' 6^ '21 



and 



si^[r + A] 



1/3 



32 



A]l/3^ 



with the discriminant 



A2 = g3 



(C9) 



(CIO) 



(Cll) 



Then the solutions can be written as 
ml = (si -f S2) - |, 

^2 = -2("'*i + '"^s) - 3 + — («i ~ ^2), 

ml = -\{s,+S2)-l-'-^{s,-S2), (C12) 

where the {777,1,7712,7773} refer to the set of masses 
{Mi,M2,M3}, but not necessarily ordered. In the one- 
loop contribution to the potential (|4.2p . only a sum over 
I enters, the order plays no role. 

We recall that the coefficients = {a, 6, c} that enter 
in the cubic eigenvalue equation \C2\ depend on the weak 
fields 77i. In the solutions (|C12p . this dependence can be 
accessed via a, si and 82- It is thus convenient to write 
(|C12I) more compactly as 



from which it follows that 
drrii 1 



daa 



X . dSr 
- Oal + ■ 

3 daa 



(C13) 



(C14) 
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Likewise, the higher derivatives are given by 



daada 



daada/jda^ 



A 



d^Sr 

daadapda^i 



(C15) 



It remains to obtain the derivatives of with respect 
to Oq,. For this purpose, it is useful to think of as a 
function of the q and r of Eqs. ((U9)) - ([UTT|l . The first 
derivatives are given by 



dsr 
daa 



dsr dq dsr dr dsr dQ" 



(C16) 



dq daa dr daa dQ^ daa ' 
where we cohectively refer to q and r as 

Q^ = {g,r}. (C17) 
In this notation, the higher derivatives can be writtes 



as 



d'^Sr 



d'^sr dQ' dQ^ dsr d'^Q' 



dttadafj 
and 

d^Sr 



dQ'^dQ*- daa dap dQ^ daadajj ' 
d^sr dQ' 9Q* 



(C18) 



daadapda^ dQ^dQ*dQ" daa dap da- 

d'^Sr 



7 



d^Q' ag* ^ d^Q" dQ' 



dQ^dQ* [daadafs da^ dapda^ duo 



d^Q' dQ* 
da-ydaa dap 



dsr 



d^Q^ 



dQ^ daadapda^ 



(C19) 



Finally, the various derivatives dQ^ /daa are given by 
dq 2 c)r I 1 ^ 



dq _ 1 
95 ~ 3' 

dc ' 



dr _ 1 

db ~ 6"' 
dr _ 1 

d^ ^ ^2' 



dsi 1 + ^ 



(C20) 



dsi ^ ^ 

dq ~ 2A^/^fTA' dr " 3 ^^^fTA' 

ds2 ^ ds2 _ 1 - i 

dq 2A^/^9^' dr 3 "V^^^ 



(C21) 



dsi -a 



da 9 A^/^/fTA 3'''^/FTaV6 9 
dsi 1 Q2 



1 



db~6A "^v^TTA 18 "V^M^' 
ds2 1 



1 

a 1-3 



56 6A'/-^/f^^ 18 "Z-^^^^' 
9^2 _ 1 1-i 
dc 6 'V^~A' 



(C22) 
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